Effects of electronic correlation on X-Ray absorption and dichroic spectra at L 2 ,3 edge 
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We present a new theoretical approach to describe X-Ray absorption and Magnetic Circular 
Dichroism spectra in the presence of e-e correlation. Our approach provides an unified picture to 
include correlations in both charged and neutral excitations, namely in direct / inversion photoe- 
mission where electrons are removed/ added, and photo absorption where electrons are promoted 
from core levels to empty states. We apply this approach to the prototypical case of L2,3 edge of 3d 
transition metals and we show that the inclusion of many body effects in the core level excitations 
is essential to reproduce, together with satellite structures in core level photoemission, the observed 
asymmetric line shapes in X-ray absorption and dichroic spectra. 
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The excitation due to X-ray absorption is a complex 
phenomenon where many body effects may play a major 
role. The independent particle picture has been shown 
to be inadequate in describing features of X-ray absorp- 
tion spectra in transition metals^ and their compounds^ 
where e-e interactions are known to be non negligible. 
To account for these failures various theoretical meth- 
ods have been developed that differ on how the e-e 
Coulomb interaction is taken into account. They range 
from atomistic approaches, where a parameterized many 
body Hamiltonian is solved via configuration interaction 
methods 2 , to solid state first-principles schemes that in- 
clude properly the structure of valence electrons^ but 
loose some of the atomic many-body effects that can be 
relevant to the physics of the process. 

In this paper we present a new approach that treats on 
the same footing the localized and itinerant character of 
electrons in a solid, exploiting on one side the first prin- 
ciple calculations of the band structure and on the other 
side the localized picture to treat many body effects as- 
sociated to electron-electron interactions. We will show 
that the on-site interaction between core and valence elec- 
trons, responsible of satellite structures in core level pho- 
toemission spectra, gives rise to asymmetric line shapes 
in XAS and XMCD spectra in agreement with experi- 
ments. Application to the prototypical case of absorption 
from the £2,3 edge in 3c? transition metals is presented. 

Within the independent-particle scheme the X-ray ab- 
sorption is described as the addition of one electron-hole 
pair to a non-interacting Fermi sea. In the language of 
many-body theory, this corresponds to approximate the 
two-particle polarization propagator to zcroth order. It 
is possible to improve this approach by substituting the 
bare particle and hole propagators with the dressed ones. 
Except for vertex corrections this would be the exact 
two-particle propagator. This picture is very physical in 
the sense that the creation of the e-h pair due to X-ray 
absorption can be visualized as the removal of one core 
electron in the presence of the valence continuum and its 
addition to the conduction band, a sort of superposition 
of photoemission and inverse photoemission spectra, plus 



possible electron-hole interaction. 

Even neglecting the excitonic effects associated to the 
electron-hole interaction we are left with a very hard task, 
namely the calculation of hole and particle propagators 
of an interacting many particle system. This is even more 
difficult here since we are dealing with two processes that 
from the point of view of e-e interactions are very differ- 
ent: on one side the removal of one electron from an 
inner core state with atomistic interactions resulting in 
multiplet structures that in most cases survive the band 
formation; on the other side the addition of one electron 
in an itinerant conduction state where local short range e- 
e interactions coexist with hopping from site to site. An 
approach that can treat these two apparently opposite 
situations in an unique way is based on the generalized 
Hubbard model where the different strength of localiza- 
tion can be taken into account in terms of different U and 
J parameters and different band widths (finite and zero 
band width for valence states and non-dispersive inner 
core energy levels respectively). In the present case the 
Hubbard Hamiltonian can be usefully partitioned in core 
(c) and valence (i>) contributions 

H = H V +H C . (1) 

The first term describes valence states in terms of single 
particle band eigenvalues (ekna) and of on-site Coulomb 
(U vv ) and exchange (J vv ) interaction 
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Here i is the site index. Similarly 

H c = H cc + H cv 

where 

iij 



2 



7: 22( U ov - JcvSo 



\s-in* r^Ti 
' J^k'a'^k 



(2) 



kk'p 00' 
a ka a k+pa l 



k' a' 



-pa' ' 



Notice that -ff c includes both core-core (U cc ) and core- 
valence (U CV ,J CV ) e-e interaction. In the present case we 
are interested in 1^2,3 edge excitations, therefore we will 
consider Coulomb and exchange integrals involving 2p 
and 3d orbitals. 

The orbital character of the valence band conies into 
our scheme through the d orbital coefficients C™(fc) of eq. 
[5]obtaincd in our case by ab-initio band structure calcula- 
tion. Both core and valence states have been calculated 
within the full-potential, linearized plane-wave method 
(FLAPW) as implemented in the Wien2k code£. Spin- 
orbit coupling has been included in order to reproduce 
a non-zero total orbital momentum that is one of the 
quantities that can be extracted from dichroic spectra. 
For core states the relativistic effects arc accounted for 
by solving the Dirac equation while for valence states the 
spin-orbit coupling is treated by perturbation theory^. 

The ab-initio band structure results could also be used 
to estimate the U and J parameters: the bare on-site 
valence - valence and core - valence coulomb and ex- 
change interactions can be explicitly calculated using the 
core and valence wavefunctions and properly scaled to 
account for screening effects. This procedure has been 
applied to the evaluation of core-core integrals in tran- 
sition metals oxides assuming a static effective dielectric 
constant^. Various ab-initio estimate of on-site valence- 
valence Coulomb and exchange integrals have been re- 
cently proposed for 3d transition metals, either within a 
constrained density functional approach^ or as dynami- 
caly screened on-site integrals^. The results depend on 
the procedure adopted and the range of variation in the 
parameters (1 to 5 eV) remains the same reported previ- 
ously in the literature^. In the present calculation U and 
J have been used as adjustable parameters, with values 
for valence-valence within the above mentioned range, 
tuned to optimize the agreement between theory and ex- 
periment. Their values are reported in table HI These val- 
ues for the valence-valence interactions give an accurate 
description of the quasiparticle band dispersion probed 
by angular resolved photoemission^—. The same is 
true of core-core interactions that have been fixed to re- 
produce the satellite position in core level photoemission 
spectra (see next figure on core level photoemission). 

When dealing with charged excitations where one elec- 
tron is removed from a core level or added to an empty 
valence state one can solve independently H c or H v to 
obtain hole and particle spectral functions respectively. 
This is done here using the 3-body scattering (3BS) ap- 
proach that has been implemented to treat both core and 
valence spectra£^~— : a charged excitation is described 
in terms of scattering between the single particle state 
with one removed / added particle and the excited config- 
urations containing one extra e-h pair. 

Energy and spin dependent complex self-energies 



TABLE I: Coulomb and exchange integrals (in eV) describing 
the interaction between valence (3d) and core (2p) states. 
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Fe 


1.5 


0.9 


0.8 


0.1 


Co 


2.1 


0.9 


1.8 


0.2 


Ni 


2.0 


0.9 


2.5 


0.2 



YfJ v (lu) are calculated and from them spectral functions 



for core and valence states given by 
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Please notice that for valence states the spectral function 
depends on band index and k- vector: A^(lo) = Akna(w) 
the k- and band index dependence coming from e v = e^ n 
and possibly from the self-energy. 

Let us start from the core hole excitation. Within the 
present approach the hole self energy turns out to b o 14 i 16 
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■[l + (U cv -J cv )]A cd (oj-e)de. (4) 



Here n da (e) is the d contribution to the valence density 
of states of spin a. The hole-hole scattering associated 
to the on-site interactions is described by T-matrices T^l 
and Tfi% for scattering between antiparallel and parallel 
spin holes respectively; similarly A a {A a ) includes an- 
tiparallel (parallel) electron-hole scattering^. 

Fig. [1] shows the calculated spectral functions for the 
2p core level of Fe, Co and Ni obtained assuming the U 
and J values of table U with binding energies adjusted 
to the experimental values^. Since our calculated self- 
energy is responsible of an intrinsic broadening identi- 
cal at the two edges, to reproduce the observed different 
lifetimes at the two edges we have added extra broad- 
ening S = Q.5eV and , S = l.heV for 2p 3 /2 and 2p 2 /2 
respectively. The spectral functions reproduce the essen- 
tial features of the core photoemission data for transi- 
tion metals where a main peak is always followed by a 
satellite structure at higher binding energies. This struc- 
ture is the fingerprint of e-e correlation being commonly 
attributed to the admixture of different single particle 
configurations: not just the ground state with one core 
electron missing but also the shake-up states where elec- 
trons are promoted to higher energy levels. These physi- 
cal processes are contained in our many body description 
based, as previously mentioned, on an extension to the 
solid state of the configuration interaction scheme. 
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FIG. 1: 2P core level photoemission for Fe, Co and Ni inclu- 
sive of correlation effects (black line) compared with experi- 
mental ones^£ (grey area). 



We turn now to valence states. The self-energy correc- 
tion to a valence band level tkna is slightly more compli- 
cated with respect to the core one reflecting the complex 
valence band structure where different orbitals are in- 
extricably hybridized (for a detailed description of the 
valence self-energy see reference^). 

The effect of correlation on empty states is shown in 
figure [2] for iron where the electron spectral functions are 
shown along high symmetry directions of the Brillouin 
zone. Notice that only some of the empty bands (the 
minority spin ones of d orbital character) are affected 
by self-energy corrections, mainly in terms of life time 
broadening. The overall effect of correlation on empty 
states is however very small, and even smaller for Cobalt 
and Nickel where d states are almost fully occupied. 

The response of interacting electrons to the creation of 
an e-h pair due to X-ray absorption is connected to these 
hole and particle spectral functions: within the single 
particle approach the absorption cross section is just the 
joint density of core and valence states modulated by the 
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FIG. 2: (color online) Iron empty states: interacting spectral 
functions are shown as color map and compared with single 
particle eigenstates reported as red (spin down) and blu (spin 
up) arrows. In this energy region self-energy corrections for 
majority spin states are negligible. 



dipole matrix elements 

M ±H (x^l^J^/mx^H (5) 
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One effect of many body correlation is to replace this 
joint density of states with the convolution of hole and 
particle spectral functions. Formally this corresponds to 
calculate a two-particle polarization propagator Xcv(w) 
where the bare particle and hole propagators are substi- 
tuted with the dressed ones 



^(W) OC \ D tkJ 2 J2 Im Xckna(uj) 
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Imxckna(u) = A c a (Cl)A kna (Cl + uj)dCl 



(6) 



with 



D ckn appearing both in cq. ((SJ) and eq. © is the 
matrix element of the electric dipole moment over the 
core and valence states 
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with the polarization vector (left/right) of the incident 
photon, e the electron charge, $J m . and the spin- 
orbit split core and valence eigenstates respectively. The 
XMCD signal is obtained as the difference /z + (w)— /j,~(lj). 

Figure [3] shows the calculated absorption (a) and 
dichroic (b) spectra for Fe, Co and Ni calculated both 
with and without self-energy corrections, compared with 
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experimental dat a 18 i 19 . We notice that in most cases 
the inclusion of many-body effects improves the agree- 
ment with experimental XAS spectra, making the lines 
more asymmetric. What makes line shapes asymmet- 
ric is the presence of shake-up satellite structures in the 
hole-spectral functions. These structures are essential to 
recover asymmetric line shapes also in XMCD spectra. 
Since in our calculation the Li and L3 lines have the 
same line shape by construction, the only difference be- 
ing their width, it is not surprising that the asymmetry 
in the dichroic lines is the same at the two edges. In the 
case of nickel where both peaks in the measured dichroic 
spectrum present the same pronounced asymmetry, the 
inclusion of e-e correlation improves the agreement be- 
tween theory and experiment. For iron and cobalt on 
the contrary we are able to reproduce rather well the 
dichroic line shape at the L2 edge but less accurately 
the I/3 one which turns out to be experimentally rather 
symmetric and more similar to the uncorrelated result. 
The asymmetry of the L2 edges increases from Fe to Ni, 
in agreement with the corresponding trend in the core- 
valence electronic correlation described by U p d (see Tab. 
!)• 

As we mentioned at the beginning what is still miss- 
ing in this description is the e-h attraction responsible of 
excitonic effects. These effects can be included by consid- 
ering the two-particle eigenvalue problem assuming the 
excited states of the N-particle interacting system to be 
a superposition of single particle states with one core 
hole and an electron above the Fermi level. This is the 
so-called Tamm-Dancoff approximation (TDA)2£ and is 
at the core also of the 3BS scheme applied here to de- 
scribe the charged excitations - just one e-h pair added 
in this case to the Nil -particle interacting system. 
TDA has been used to solve the Bethc-Salpeter equa- 
tion (BSE) and treat excitonic effects in valence- valence 
transitions^ and more recently extended to X-ray ab- 
sorption from core states^^. In both cases the dressed 
one-particle propagators have been obtained within the 
GW approximation^. In the present case we are inter- 
ested in a regime of high e-e correlation where the a per- 
turbativc approach like GW is supposed to be less ap- 
propriate being unable to reproduce satellite structures 
that dominate the core level spectra. A non-pcrturbative 
method such as the 3BS approach would be more appro- 
priate and we are presently testing the extension of our 
theory in this direction. 

A further difference between the present approach and 
the ab-initio methods based on cither GW+BSFj^2^ or 
time-dependent local density approximation^ is the treat- 
ment of the e-e interaction: in these schemes the elec- 
trons interact through the full long range, dynamically 
screened Coulomb potential while in our approach only 
the short-range part of the Coulomb interaction is consid- 
ered, in the spirit of the Hubbard model. The implicit as- 
sumption, common to most approaches used to describe 
e-e correlation in strongly correlated materials, is to at- 
tribute many body effects to just the on-site interaction, 
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FIG. 3: Absorption (left panels) and dichroic (right panels) 
spectra of Fe, Co, Ni. Theoretical spectra, with (QP) and 
without (SP) self-energy corrections are compared with ex- 



perimental ones 



represented by gray areas. 



the long range part being included at the mean-field level 
in the band calculation. This assumption seems particu- 
larly reasonable in the case of the strongly localized core 
states but in practice leads to the introduction of ad- 
justable parameters. This is has been up to now a short- 
coming of most the approaches based on the Hubbard 
model, with the advantage to allow for solutions beyond 
the perturbation limit. 

In summary, we have shown that many-body effects 
due to e-e interactions can be included in the description 
of X-ray absorption and dichroic spectra by replacing the 
joint density of states appearing in the absorption cross- 
section with a convolution of (core) hole and (valence) 
electron spectral functions. Whenever e-e correlations 
modify significantly hole and electron spectral functions 
we expect non negligible effects in the absorption spectra 
as well. In the case of 3d transition metals these modifica- 
tions are essentially associated to the shake-up structures 
that appear in the core-level photoemission spectra as a 
residue of atomic multiplets. These structures, related to 
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the admixture of different single particle configurations, 
can be reproduced only in a many body approach and 
are responsible of the asymmetric line shapes of XAS 
and XMCD structrues at L2 ; 3 edges in Fe, Co and Ni. 
In other systems where local e-e correlations affect sig- 
nificantly the empty part of the valence band we expect 
even more drastic changes: this should be the case for 
instance of highly correlated materials where on site e-e 
repulsion is responsible of the creation of extra empty 
(Hubbard) bands^. 

Finally we want to stress that our approach has a 
purely experimental counterpart: it would be interest- 
ing to check on an experimental basis wether - and to 
what extent - the measured absorption spectra can be re- 



produced as a straightforward convolution of direct and 
inverse photoemission data, in this way confirming the 
direct link between neutral excitations, due to X-ray ab- 
sorption, and charged excitations corresponding to hole 
and electron addition, link that is clear theoretically but 
that has not been demonstrated experimentally so far. 
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